Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 11/12/01 



SPECTRAL ENERGY DISTRIBUTION IN SUPERCRITICAL DISK ACCRETION FLOWS 
THROUGH PHOTON-TRAPPING EFFECTS 

K. Ohsuga, S. Mineshige, and K. Watarai 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
Draft version February 2, 2008 

ABSTRACT 

We investigate the spectral energy distribution (SED) of the supercritical disk accretion flows around 
black holes by solving the multi-frequency zeroth moment equation of the radiation transfer equation 
under the flux-limited diffusion approximation as well as the energy equation of gas. Special attention 
is paid to the photon-trapping effects, the effects that photons are trapped within accretion flow and 
are swallowed by a black hole, with little being radiated away. It is found that when flow luminosity is 
below about twice the Eddington luminosity [L <2Le) photon trapping is ineffective and the SED shifts 
to the higher-energy side as L increases. The peak frequency at which the SED reaches its peak becomes 
at maximum three times higher than that given by the standard-disk model, which might resolve the 
so-called 'too hot accretion disk problem.' When L>2Le, in contrast the peak frequency of the SED 
tends to decrease with increase of mass-accretion rate. This is due to enhanced photon-trapping. Since 
high-energy photons are generated near the equatorial plane, they can be more effectively trapped in flow 
than low-energy ones and, hence, the high energy part of radiation is suppressed. Finally, the observed 
time variation of IC 342 SI, which is an ultraluminous X-ray source, in the X-ray HR diagram can be 
explained by the modulation of the mass accretion rate. In comparison with the observational data, we 
estimate the mass of black hole in IC 342 SI is ~ IOOMq. 

Subject headings: accretion: accretion disks — black hole physics — radiative transfer 



1. INTRODUCTION 

The standard-disk model proposed by Shakura & Sun- 
yaev (1973) has had some success in describing optically 
thick accretion flow structure in luminous objects, such 
as X-ray binary sources in the soft state. However, the 
standard-disk model is constructed under number of ap- 
proximations and its treatment is only justified when the 
mass-accretion rate, M, is less than the critical mass- 
accretion rate, M crrt , where the critical mass-accretion rate 
is given by Merit = Le/c 2 with the Eddington luminosity, 
Le = 4^cGMm p /aT, where c, M, m p , and ctt are the 
light velocity, the black-hole mass, the proton mass, and 
the Thomson scattering cross-section, respectively. Since 
the mass- accretion rate is determined by some mechanism 
which is independent of the nature of the central star it- 
self (Alme & Wilson 1976; Spruit & Ritter 1983; Norman 
k Scoville 1988; King et al. 1997; Umemura, Fukue, & 
Mineshige 1998; Ohsuga et al. 1999), there are no rea- 
sons why M should always be below M cr ; t . That is, the 
supercritical accretion flow may naturally arise. Indeed, 
several possible sites of super- or near-critical accretion 
flow have been suggested; e.g., micro-quasars (Yamaoka, 
Ueda, & Inoue 2001), ultraluminous X-ray sources (ULXs) 
(Okada et al. 1998; Colbert & Mushotzky 1999; Mizuno 
et al. 1999; Makishima et al. 2000; Watarai et al. 2000; 
Revnivtsev et al. 2002), narrow-line Seyfert 1 galaxies 
(Boiler 2000; Mineshige et al. 2000), quasars (Collin et al. 
2002), and Gamma-ray bursts (Narayan, Piran, & Kumar 
2001; Kohri & Mineshige 2002). 

When the mass-accretion rate exceeds the critical value, 
photon trapping plays an important role within the trap- 



ping radius, 

rtrap ~ m r s , (1) 

where rh is the mass-accretion rate normalized by the crit- 
ical mass- accretion rate, H is the disk-half thickness, r is 
the radius, and r g is the Schwarzschild radius (Ohsuga et 
al. 2002, hereafter referred to as Paper I, see also Begel- 
man 1978 for the spherical case). Thus, the radiation en- 
ergy, which is generated by viscosity in the vicinity of the 
black hole, is trapped in the flow and advected inward 
without being immediately radiated away. Since the ad- 
vective energy transport was firstly considered in the slim- 
disk model by Abramowicz et al. (1988), in the context 
of the accretion-disk models this model has been believed 
to be the solution of the supercritical disk accretion flow 
(Kato, Fukue, & Mineshige 1998). However, it is shown in 
Paper I that the photon-trapping effects are not accurately 
treated in the slim-disk approach. The reason for this is 
that the time delay between the energy generation deep 
inside the disk and the energy release at the surface is not 
taken into account in the slim-disk formulation, although 
the time delay is not negligible and exceeds the accretion 
timescale within the trapping radius. 

In the light of the recent discovery of the supercriti- 
cally accreting objects, it is of great importance to study 
the spectral energy distribution (SED) of the supercrit- 
ical disk accretion flow. But, little is known about the 
SED of the flow suffering the photon-trapping effects ex- 
cept for the case of the slim-disk model. Szuszkicwicz, 
Malkan, & Abramowicz (1996) and Watarai et al. (2000) 
have calculated the SED of the slim disk by superposi- 
tion of blackbody/modified-blackbody spectra. Wang et 
al. (1999) have solved the radiation transfer equations in 
the vertical direction, coupling with the slim-disk formu- 
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lation. As we noted above, however, the slim-disk model 
has not fully considered the photon-trapping effects, so we 
need further studies. 

In paper I, we investigated the photon-trapping effects in 
the disk accretion flow, but only solved the monochromatic 
radiation transfer equations inside the accretion flow. In 
this work, we investigate the SED of the optically-thick 
supercritical accretion flows by solving multi-frequency ra- 
diation transfer equations. In §2, we give our model and 
the basic equations. The resultant SEDs are presented in 
§3. Comparing our results with the observational data, 
we discuss the role of the photon trapping in the actual 
objects. Finally, §4 and §5 are devoted to discussion and 
conclusions, respectively. 

2. OUR MODEL AND BASIC EQUATIONS 

In this section, we numerically solve the multi-frequency 
radiation-transfer equations, as well as the energy equa- 
tions of gas and radiation, by taking into account the 
photon-trapping effects. Since we are concerned with the 
radiation processes, especially the photon-trapping effects 
themselves, we employ a simple model for the accretion 
flow structure (see Paper I). 

Here, the cylindrical coordinates, (r,(p,z), are used. 
We consider that the accretion disk is axisymmetric and 
steady in the Eulerian description; d/dtp — d/dt — 0. We 
employ the Gaussian density profile in the vertical direc- 
tion, 

T / v. \ 21 

p(r,z) = poO)exp 



The vertical component is prescribed as 



H 



(2) 



where po is the gas density on the equatorial plane and we 
assume the ratio of the disk-half thickness to the radius, 
H/r, to be a constant in radius, depending on m; 

with £ioo being a dimensionless parameter representing 
the aspect ratio, H/r, at m = 100. Such a linear rela- 
tion between H/r and m exactly holds in the radiation 
pressure-dominated part of the standard-disk model [see, 
e.g., equation (3.61) of Kato, Fukue, & Mineshige (1998)], 
and also roughly holds in the slim-disk model as long as 
rh is less than a few hundreds (Watarai et al. 2000). The 
density on the equatorial plane, po, is related to the surface 
density as, 



](r) = 2p a (r) f 



exp 



H 



dz = 2^p (r)H(r), 



(4) 

where the surface density, E, is calculated by the continu- 
ity equation for a given accretion rate, M, 



M = -2TTrv r Y„ 



(5) 



with v r being the radial component of the velocity. 

Since the radial component of the velocity, v r , is related 
to the viscosity parameter a through v r ~ — a (H/r) c s 
with c s being the sound velocity (see, e.g., Frank, King, 
& Raine 1985; Kato, Fukue, & Mineshige 1998), and H is 
given by the equation (3), v r is finally expressed in terms 
of the free-fall velocity as 
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(6) 



v z = —v r 
r 



(7) 



i.e., we assume a convergence flow. 

Through equations (3), (5), and (6), E and p are re- 
lated to £ioo and m in our model; E oc e^m -1 and 



p cx e 10 o r " • The photon-trapping radius is rewritten 
by using equation (3) as, 

rtrap ~ 10~ 2 £ioo™ 2 r g . (8) 

Here, it is noted that the photon-trapping radius not only 
indicates the critical radius for the photon trapping but 
also expresses the significance of the photon-trapping ef- 
fects, since this relation is derived from the comparison 
between the accretion timescale and the radiative diffu- 
sion timescale. Thus, the photon-trapping effects are more 
pronounced in the case of larger £i o and rh. 

We suppose the structure of the accretion disk to be 
locally plane parallel; that is, the disk is regarded as con- 
sisting of accreting ring elements (see discussion in the 
end of this section). Each ring is composed of N lay- 
ers in the z-direction, and the thickness of each layer is 
Az — 2Hq/N (see Figure 2 of Paper I). Here, Hq defines 
the upper boundary of the disk in the numerical simula- 
tions in such a way that p(r, Hq) — O.lpo(r). We solve 
the time-dependent energy fields of gas and radiation in 
the ring element during the course of accretion motion un- 
til the element reaches the inner edge of the disk. Both of 
the viscous heating and comprcssional heating though con- 
verging inflow are considered. Since the flow assumed to 
be steady (d/dt = 0), time coordinate can be transformed 
to the spatial coordinates; i.e., D/Dt — v r d/dr + v z d/dz 
with v r and v z being given by equations (6) and (7). We 
thus express the multi-frequency radiative flux at the disk 
surface as a function of radius and then calculate the re- 
sultant SED and luminosity of the disk. 

In the plane-parallel approximation, the radial and az- 
imuthal components of the radiative flux are null. Also, 
non-diagonal components of the radiation stress tensor are 
null. Hence, using equations (6) and (7), we write the en- 
ergy equations of radiation and gas as 
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and 
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= - ~ — Pgas - 4:irKpB +cn E E + q vis , (10) 

A V 

respectively, where E v and E(= J E v dv) are the radi- 
ation energy density and its frequency-integrated form, 
respectively, F* is the radiative flux in the z-direction, 
P™ is the diagonal (n-th) components of the radiation 
stress tensor, B v and B = ctT^/tt are the Planck func- 
tion and its frequency-integrated form, respectively, with 
T gas being the gas temperature and a being the Stefan- 

Boltzmann constant, k v = 3.7 x 10 8 T ga s / ' 2 (p/mp) 2 i/ _3 {l - 
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exp(—ft,^/fcT gas )} is the free-free absorption coefficient 
(Rybicki & Lightman 1979), with h being the Planck con- 
stant and k being the Boltzmann constant, e is the internal 
energy density, p gas is the gas pressure, Kp and ke are the 
Planck-mean and energy-mean opacities, respectively, 



Kp 



jjj J K V B V 



dv, 



— I K v E v dv 1 



(11) 
(12) 



and qvis is the viscous heating rate per unit volume (Miha- 
las & Klein 1982; Mihalas & Mihalas 1984; Fukue, Kato, 
& Matsumoto 1985; Stone, Mihalas, & Norman 1992). 

To close the set of equations, we apply the multi- 
frequency flux- limited diffusion approximation (FLD), for 
the radiative flux and stress tensor, 



F z = -- 



Xv 

1 



dz 



and 



pz 



with Xv being the extinction coefficient, 

pcr T 

Xv = h K v , 

m p 

where X v , R v , and f u are defined as 
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(13) 
(14) 

(15) 

(16) 

(17) 
(18) 



and 



f„ = K + \ 2 v R 2 v . (19) 

These relations correspond to the multi-frequency form of 
the monochromatic FLD represented by Turner & Stone 
(2001). This approximation holds both in the optically 
thick and thin regimes. In the optically thick limit (R v — > 
0), we find X v — ► 1/3 and f„ — > 1/3. In the optically thin 
limit of i?„ — ► oo, on the other hand, we have \F£\ = cE u , 
prr _ _ q ; and pzz _ p^ These give correct re- 
lations in the optically thick diffusion limit and optically 
thin streaming limit, respectively. Thus, we can rewrite 
the radiation energy equation (9) as 
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dr 
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E, 



'dz \Xv dz 
+Akk v B v -ck„E v . (20) 
These nonlinear equations (10) and (20) are integrated 
iteratively by the Ncwton-Raphson method with the 
Gauss- Jordan elimination for a matrix inversion, by be- 
ing coupled with the equation of state, 

2 

Pgas = gC (21) 

We employ the simple model for viscous heating rate, 
whose distribution in the z-direction is proportional to the 
total pressure; that is, 

E(r, z) + 2e(r, z)/3 



q v is(r, z) = Q vis (r) 



C{E(r,z) + 2e(r,z)/3}dz' 



(22) 



where Q v is is the half vertically-integrated viscous heating 
rate per unit surface, 

1/21 



8tt 



1 - 



(23) 



(Shakura & Sunyaev 1973; Lynden-Bell & Pringle 1974), 
with f^K being the Keplerian angular speed. 

Throughout the present study, a and N are set to be 0.1 
and 100, respectively. We consider the range of frequency, 
logi^ = 16 — 20, with Alogf = 0.08. Since the radia- 
tive flux from the vicinity of the black hole gives negligible 
contribution to the total SED due to the large redshift in 
the face-on case, we neglect the emission from the regions 
inside, r < 3r g [see, e.g., equation (3.159) of Kato, Fukue, 
& Mineshige (1998)]. The supercritical disk becomes rel- 
atively geometrically thick through the strong radiation 
pressure and is obscured by the outer part of the disk in 
the edge-on case so that it would be observed as objects of 
small inclination angle. The black-hole mass is fixed to be 
1OM except in §3.4, where we adopt 5OM and 100M Q 
to compare our results with the observations of ULXs. 

In our simulations, we assumed the disk to be locally 
plane parallel. Hence, the radiation energy diffuses to- 
wards only the z-direction in the frame moving with gas; 
i.e., the radial and azimuthal components of the radiative 
flux arc null. However, if the radial gradient of the ra- 
diation energy density is comparable to the vertical one, 
\dE/dz\ ~ \dE/dr\, the diffusion of the radiation energy 
towards the r-direction is not negligible. Such a situation 
would hold in the very vicinity of the inner edge of the 
disk, where the radiation energy effectively falls onto the 
black hole because of F r < 0. Thus, the photon trapping 
would be enhanced at r <3r g . In contrast, the radiative 
flux would become positive, F r > 0, at r > 3r g , and may 
not be negligible if the aspect ratio, H/r, is of the order 
of unity. In this case, the photon-trapping effects would 
be reduced. Such a reduction of the photon trapping is 
roughly estimated as follows. Since the radiative diffusion 
timescale is given by c/t, the ratio of the radial diffusion 
velocity and vertical diffusion velocity is roughly obtained 
to be Vdifl\r/vdiff,z ~ t z /t t ~ H/r. Here, r z ~ paTH/m p 
and r r <~ p<TTr/m p are the optical depth for the Thomson 
cross-section in vertical and radial directions, respectively. 
By taking into account the radiative diffusion in the ra- 
dial direction, the timescale for the radiation energy to 
be swallowed by the black hole is r swa ~ r/(v r — Wdiff.r)- 
By comparison between this timescale and the diffusion 
timescale in vertical direction, H/(c/r z ), we find that the 
photon-trapping radius is modified as 



trap 



^trap 




(24) 



where equations (1), (5), r z = or^/mp, M cr ; t = Le/c 2 , 
and Le = 4ircGMm p / '<jt are used. Since the luminosity 
depends on the trapping radius as L <~ f 2irrQ v i s dr oc 

r tra P (Paper I), the luminosity would increase by a factor 
of 1 + (H/r) 2 on account of the radiative diffusion in the 
radial direction. For eioo = 0.5, our formulation under- 
estimate luminosity by (H/r) 2 <~ 0.25(m/100) 2 . Dopplcr 
beaming tends to focus radiation field to be more confined 
within the equatorial plane, thereby enhancing the photon 
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trapping effects. The details would be elucidated by solv- 
ing two-dimensional radiation transfer equations as future 
work. 

3. RESULTS 

3.1. Luminosity 

We first plot the luminosity against the mass- accretion 
rate in Figure 1 , where the luminosity and mass-accretion 
rate are normalized by the Eddington luminosity and the 
critical mass-accretion rate, respectively. The luminos- 
ity should increase along the dashed line, if the energy- 
conversion efficiency, r\ = Mc 2 /L, is constant. How- 
ever, it is found that the energy-conversion efficiency 
decreases with increase in the mass-accretion rate at 
L >2Le due to the photon trapping. We also see that 
the larger eioo is, the lower becomes the luminosity at 
high mass-accretion rates, rh ^> 10. This can be un- 
derstood, since the photon-trapping effects are more con- 
spicuous for larger eioo [see equation (8)]. We have con- 
firmed the results of Paper I, whereby energy-conversion 
efficiency decreases with increase of the mass-accretion 
rate, by the multi-frequency radiation transfer simulations. 
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Fig. 1. — Total disk luminosity as a function of the mass-accretion 
rate, rh = Mc 2 / Lg. The circles, squares, and triangles represent the 
cases with eioo = 0.5, 0.17, and 0.1, respectively, where eioo is a 
dimcnsionless parameter denned as eioo = (H/r)/ (rh/ 100). The 
dashed line shows the luminosity for the cases that the energy- 
conversion efficiency, r) = Mc 2 / L, is kept constant. 

3.2. Effective Temperature 

The effective temperature profile is shown in Figure 
2. Here, the effective temperature is defined as T e g- = 
(F/cr) 1 / 4 , where F = J F v dv. The profile becomes a 
bit flatter in the vicinity of the black hole when m >30 

due to the photon-trapping effects. In particular, the 
slope becomes somewhat flatter than that predicted by 
the slim-disk model, r^ 1 / 2 , in the case of rh = 100. When 
rh = 10, in contrast, the profile is consistent with the 
standard-disk model, in which T c s oc r~ 3 / 4 , except at 



r/r g <7, where boundary effects appear, since the photon- 
trapping effects are negligible. Here, it should be stressed 
that the flatter effective temperature means a reduction 
of the radiative flux at small radii, however, the lo- 
cal emergent spectrum is no longer a simple multi-color 
blackbody spectrum (discussed in the next subsection). 
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Fig. 2. — The effective temperature profiles for the case with 
£ioo = 0.5 for various rh; rh = 10, 30, and 100. Through the 
photon-trapping effects, the profile becomes flatter in the vicinity 
of the black hole as m increases beyond m ~ 30. In particular, 
the profile for the case with rh = 100 is somewhat flatter than the 
prediction of the slim-disk model; i.e., T e g oc r -1 / 2 . 

3.3. Spectral Energy Distribution 

We show the SED of the accretion flow with m = 10 
for £ioo = 0.5, 0.17, and 0.1 in Figure 3. For such a low 
rh, the photon-trapping effects do not appear yet, but the 
calculated SEDs grossly deviate from the SED predicted 
by the standard disk with the same rh, which is indicated 
by the dotted curve. It is found that the SED is shifted 
towards the higher-energy side as eioo increases. This is 
explained by the transmission of high-energy photons gen- 
erated near the equatorial plane, since for given rh both of 
E and p decrease as eioo increases, leading to a decrease of 
k v , especially at high v (note k v cx p 2 and k„ decreases as 
v increases). Then, the high-energy photons arising from 
deep inside the disk can pass through the disk body when 
eioo is large. In the case of eioo = 0.5 and 0.17, the fre- 
quency of the SED peak intensity (hereafter referred to as 
the peak frequency) is by about three times higher than 
that of the standard disk. This result might resolve the 
so-called 'too hot accretion disk problem' (Makishima et 
al. 2000). That is, the black-hole masses of some ULXs 
are inferred to be at least ~ IOOMq, since their luminosi- 
ties reach ~ 10 40 erg s _1 , while the source luminosities can 
not exceed the Eddington luminosity. Then, the standard- 
disk model predicts that the innermost disk temperatures 
of these objects should be kT in < 1 keV for M >100M Q , 

however, the observed temperatures are much higher (1.1- 
1.8 keV). This is the 'too hot accretion disk problem.' 
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Fig. 3. — The spectral energy distribution (SED) of the accre- 
tion disk for rh = 10 and for eioo = 0.5, 0.17, and 0.1. The dotted 
curve shows the SED predicted by the standard-disk model. The 
SED is shifted towards the higher-energy side with the increase of 
eioo, since for large £ioo,S, p, and, hence, k v are smaller so that 
the high-energy photons arising from deep inside the disk can easily 
pass through the disk body. 

In our simulations, the spectral hardening of the emer- 
gent SED is caused by the transmission of the high-energy 
photons arising from deep inside the disk. The Comp- 
tonization is not taken into account in our simulations al- 
though it provides another important cause of the hard- 
ening. The hardening due to the Comptonization was 
demonstrated by Czerny & Elvis (1987) and Ross, Fabian, 
& Mineshige (1992) [see, however, Laor & Netzer (1989) 
for negative results]. Shimura & Takahara (1993) solved 
multi-frequency radiation transfer equations coupled with 
the hydrostatic balance and Comptonization, finding the 
spectral hardening factor of n <~ 1.7. 

When eioo is small, say eioo = 0.1, S and 
large, and, therefore, photons tend to be absorbed and re- 
emitted before reaching the disk surface. Thus, only low- 
energy photons which are generated near the disk surface 
can go out and thus the emergent SED becomes similar to 
that of the standard disk. However, since the absorption 
coefficient by the free-free absorption strongly depends on 
the frequency as a part of high-energy photons can 
directly pass through the disk, and a double peak appears 
(see the cases with eioo = 0.17 and 0.1). Here, we should 
note that the higher-energy peak might disappear, if the 
metal absorption, which can be much more effective at 
high frequency, is included (discussed later). 

The emergent SEDs in the case of rh — 100 is 
shown in Figure 4. As was mentioned above, the 
peak frequency tends to increase as eioo increases due 
to the transmission of the high-energy photons. How- 
ever, the SED in the case of eioo = 0.5 is softer than 
other cases. This is caused by the photon-trapping ef- 
fects. Since the higher-energy photons originate from 
deep inside the disk, they tend to be more effectively 
trapped in the flow than lower-energy ones. As a con- 



sequence, the peak frequency of the emergent SED even- 
tually turns to decrease as rh further increases. That 
is, the photon trapping provides the counter effect to 
the transmission of the high-energy photons. Thus, the 
peak frequency of emergent photons starts to decrease 
when i7i exceeds some critical value [see equation (8)]. 
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Fig. 4. — Same as Figure 3, but for rh = 100. The thin solid 
curve indicates the SED of the multi-color disk (MCD) model. The 
peak frequency of the SED becomes considerably low in the case of 
eioo = 0.5 due to effective photon trapping. The MCD model does 
not correctly reproduce the actual SED. 

In Figure 4, the thin solid curve is the SED of the multi- 
color disk (MCD) model for eioo = 0.5, where the SED 
of the MCD model is calculated by the superposition of 
the blackbody spectra with various effective temperature, 
B v {T eS ) with T cff = (F/a) 1 /*. As shown in this figure, the 
emergent SED is deviated from that of the MCD model. 
Therefore, not only the standard-disk model but also the 
MCD model can not successfully reproduce the actual 
SEDs in the supercritical accretion flow, and it is nec- 
essary to incorporate photon-trapping effects accurately 
when rh >100. 

In Figure 5, we plot the variations of SED with changes 
in rh for eioo = 0.17 and 0.5. Symbols filled squares and 
circles indicate the locations of the SED peak. The peak 
frequency of the SED increases as the mass-accretion rate 
increases up to rh = 30 for both cases. When rh is larger 
than 30, however, the peak frequency no longer increases, 
since the photon trapping works and tends to suppress 
the emergence of high-energy photons. In the case of 
£ioo = 0.5, in which the photon-trapping effects are more 
enhanced than the case of ei o = 0.17, the peak frequency 
rather decreases with increase of the mass-accretion rate. 
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Fig. 5. — The SED for eioo = 0.17 (upper panel) and 0.5 (lower 
panel) for different mass-accretion rates, rh = 10, 30, and 100. In 
contrast with the standard disk regime at rh < 30, in which the 
peak frequency increases with increasing rh, the peak frequency of 
the SED does not increase with an increase of mass-accretion rate 
when the photon trapping is effective (rh > 30). Rather, the peak 
frequency decreases as rh increase, when eioo = 0.5. 

To sum up, the SED of the supercritical accretion flow 
shifts towards the higher-energy side owing to the trans- 
mission of high-energy photons as long as the photon- 
trapping effects are not appreciable. When rh exceeds 
<~ 30, however, the photon-trapping effects become sub- 
stantial if £ioo^>0.17, and the peak frequency starts to 
decrease in spite of an increase of the mass-accretion rate. 

3.4. Spectral fitting 
3.4.1. X-ray H-R diagram 

In the previous subsection, we overviewed the emer- 
gent SED in the supercritical accretion flow with care- 
ful considerations on the photon-trapping effects. Here, 
we fit the calculated spectra, following Watarai et al. 
(2000). In this fitting method, a simple temperature pro- 
file, Ti n (r/ri n ) _3 / 4 with T- m and r- m being fitting parame- 
ters, is employed for the model SED (Mitsuda et al. 1984). 
We adopt the photon energy range between O.lkeV and 
5keV for the fitting, but we should keep in mind that 
the fitting results depend on the adopted energy range, 
in which the fitting is made. We plot the calculated T ln as 
a function of L in comparison with the observational data 
in the X-ray H-R diagram (see Figure 6). 

As shown in Figure 6, kT- m increases as the mass- 
accretion rate increases for the case of M = 10M Q and 
£ioo = 0.17 (filled squares). In the case of eioo = 0.5 
(filled circles), in contrast, the innermost temperature Tj n 
decreases with increase of rh when rh > 30. This is 
due to very efficient photon-trapping effects, by which 
the emergent SED shifts towards the lower-energy side, 
leading to a decrease in kT- ln . As was mentioned in 
the previous subsection, the photon-trapping effects are 
more conspicuous in the accretion flow with larger eioo- 
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Fig. 6. — The predicted and observed trajectories of black-hole 
object with varying luminosity in the X-ray H-R diagram. The filled 
circles and filled squares indicate our numerical results for eioo = 0.5 
and 0.17, respectively, with different mass-accretion rates, rh = 10, 
30, 50, and 100 from the bottom to the top. The small and big 
open squares represent the results for eioo = 0.17, but the black- 
hole mass is M = 50Mq and IOOMq, respectively. A result of 
rh = 100 is not shown for M = 50Mq, and we plot only results of 
rh = 10 and 30 for M = IOOMq. The dotted lines correspond to 
the cases of the standard-disk model with different black-hole mass, 
M = 1M , 3Mq, WMq, and 33M . According to the standard- 
disk model, the luminosity becomes the Eddington luminosity on 
the thin dashed line. The case of (M, eioo) = (10Mq,0.17) nearly 
coincides with the observational results of M81 S-6. Moreover, the 
observed variation of IC 342 SI can be reproduced by our model for 
(M, eioo) = (50 - 1OOM , 0.17). 

In this figure, we found that the case of eioo = 0.17 
(filled squares) coincides most with the observational re- 
sults of M81 S-6. The cases of eioo = 0.5 (filled circles) 
does not reproduce the observed trend. Therefore, we 
can rule out the cases of eioo = 0.5, that is, the photon- 
trapping effects should only moderately appear in the ac- 
cretion flow of ULXs. This conclusion is supported by 
the fitting to the IC 342 SI data, which is also well re- 
produced by the cases with eioo = 0.17, if the black-hole 
mass is scaled up to M = 50 — 100M Q (small and big open 
squares) . 

From the fitting to the IC 342 SI data, we find that the 
black hole in IC 342 SI is considerably massive; its mass 
is 50 — 100M©. Note that our estimation is by a factor 
of 1.7 — 3 larger than the value ~ 30M© which Watarai, 
Mizuno, & Mineshige (2001) estimated based on the slim- 
disk model. Scatters of IC 342 SI data in the H-R diagram 
is understood in term of the variation of the mass-accretion 
rate according to our model. If we adopt the standard-disk 
model, the black-hole mass should change to fit the obser- 
vation which is, of course, unreasonable. Our model is 
certainly preferable. 

3.4.2. Inner-edge radius and hardening factor 

Next, we estimate the inner-edge radius, R ln , by adopt- 
ing equation (4) in Makishima et al. (2000) to our cal- 
culated SEDs. We assume that the hardening factor is 
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k = 1.7 and the correction factor to i?; n , which is re- 
quired to remove boundary effects is set to be £ = 0.42. 
As shown in the upper panel of Figure 7, the resultant 
inner-edge radius depends on the parameter, eioo, as well 
as the mass-accretion rate, to; namely it depends on how 
important the photon-trapping effects are. The photon- 
trapping effects tend to increase the inner-edge radius 
and decrease of T- m because of suppression of emergence 
of high-energy photons. In the case of eioo = 0.5, i?; n 
drastically shifts from 40 km to 160 km as the mass- 
accretion rate increases from to = 30 to 100 due to very 
efficient photon trapping. If £i o = 0.17, conversely, 
the radius decreases with increase in the mass-accretion 
rate as long as m<50. At higher to > 50, the radius 
slightly increases due to the mild photon-trapping effects. 
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Fig. 7. — The inner-edge radius, R ln (upper panel) and the hard- 
ening factor, k (lower panel) for the different mass-accretion rates, 
fa = 10, 30, 50, and 100, from the left to the right, respectively. 

In the lower panel of Figure 7, we plot the hardening 
factor, «, by employing equation (4) in Makishima et al. 
(2000). Here, we set R ia = 3r g and £ = 0.412. It is 
clear in this figure that the hardening factor is reduced by 
the photon-trapping effects. In the case of £ioo = 0.5, ft 
decreases when to exceeds 30. If £ioo = 0.17, k slightly 
increases for to < 50, but it decreases at to > 50 when 
photon-trapping effects are substantial. 

4. DISCUSSION 

4.1. Comparison with the Slim-disk Model 

In some slim-disk approach, the emergent SED is as- 
sumed to be the superposition of blackbody spectra with 
different effective temperatures (Szuszkiewicz, Malkan, & 
Abramowicz 1996; Watarai et al. 2000). However, such a 
method is not successful in reproducing the actual SED as 
was shown in §3.3. Moreover, we found that the effective 
temperature profiles may deviate from the prediction of 
the slim-disk model, r -1 / 2 , at large to since the photon- 
trapping effects are not accurately treated in the slim-disk 
model (§3.2, see also Paper I). 



Nevertheless the slim-disk model, as well as our model 
for the case of £ioo = 0.17, can reproduce the basic trend 
of the ULXs data in the X-ray H-R diagram (Watarai, 
Mizuno, & Mincshigc 2001). However, we should stress 
that the physical reason to cause apparent decrease in 
i?in (or ratio of i/T 4 ff ) completely differ. In the slim- 
disk approach, the physical inner edge actually shifts and 
extends down to ~ r g at high m. The consequence is 
that the ULXs data can be fitted well with small R ln . In 
our model, in contrast the disk inner edge is assumed not 
change, but the disk temperature increases apparently due 
to the transmission of the high-energy photons from deep 
inside. In realistic situations, both will work, but which is 
more important is not clear at the moment. 

Finally, we indicate the important difference between 
the slim-disk model and our model. The black-hole mass 
of IC 342 SI is estimated to be ~ 3OM in the slim- 
disk approximation (Watarai, Mizuno, & Mineshige 2001), 
whereas it is considerably massive, 50 — 100M Q , in the 
present study. We need to remark on the limitations of 
the present model, however, since our model is essentially 
one-dimensional model and the flow structure is not prop- 
erly solved. We need further studies (see below). 

4.2. Future Work 

Throughout the present study, we treated only the free- 
free absorption/emission and the Thomson scattering as 
physical interaction between gas and radiation. However, 
the Comptonization may play an important role on the 
emergent SED (Czerny & Elvis 1987; Ross, Fabian, & Mi- 
neshige 1992; Shimura & Takahara 1993). Then, the emer- 
gent spectra get even harder. If the corona exists above 
the disk, farther, the high-energy power-low photons will 
be generated by the inverse Compton scattering (Haardt 
& Maraschi 1991; Liu, Mineshige, & Ohsuga 2002). The 
emergent SED would be improved by simulations of ra- 
diation transfer which incorporates Comptonization, e.g., 
Monte Carlo simulations (Pozdnyakov, Sobol', & Sunyaev 
1977). 

The metal absorption may be effective at high frequency. 
Then, the high-energy photons generated deep inside the 
disk may be absorbed at a certain depth. If this is the 
case, the higher-energy peak may disappear and the peak 
of the emergent SED shift to lower-energy side. We ex- 
pect that the emission lines by metals should be observed 
in the emergent SED. If the cross-section for the metal 
absorption is comparable to or more than that of the 
Thomson scattering, the photon trapping is drastically 
enhanced and the energy-conversion efficiency is more re- 
duced. Then, such accretion flows will be observed as rel- 
atively dim and metal rich objects. 

In this study, we adopted the plane-parallel approxima- 
tion, which introduces error of ~ (H/r) 2 <C 1. As fu- 
ture work, we need to solve full two-dimensional radiation 
transfer equations. 

The convection which seems to occur when the disk is 
radiation pressure dominated (Shakura, Sunyaev, & Zil- 
itinkevich 1978; Agol et al. 2001) would also reduce the 
photon-trapping effects. Since the high-energy photons 
trapped in the flow can be carried towards the disk sur- 
face by the convection, the high-energy photons emerge 
and the peak of the SED would be higher than that of 
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the expectation in our present study (see discussion in Mi- 
neshige et al. 2000). 

The photon-trapping effects would be much reduced, 
if the gas density distribution is strongly patchy in the 
disk, since the photons can pass through the disk along 
the low-density regions. Begelman (2002) suggested that 
by such a effect the disk luminosity can exceed the Ed- 
dington luminosity. The details are not clear before multi- 
dimensional radiation-hydrodynamical (RHD) simulations 
are performed. Two-dimensional RHD simulations have 
been initiated by Eggum, Coroniti, & Katz (1987, 1988), 
and followed by Okuda, Fujita, & Sakashita (1997), Fu- 
jita & Okuda (1998), Kley & Lin (1999), Okuda & Fujita 
(2000), and Okuda (2002), though the radiation fields are 
treated in monochromatic form. 

Finally, wc employ the FLD approximation in the 
present work. Though this formulation well agrees with 
the exact solution both in the optically thick and thin 
regimes, it would be better to solve the radiation transfer 
equations without using the approximation. Especially, 
the FLD approximation might be invalid when the patchy 
structure forms inside the accretion flow, since the FLD 
approximation does not always give the exact solution at 
the regions of moderate optical thickness. 

5. CONCLUSIONS 

By solving the multi-frequency radiation transfer and 
the energy equations of gas as well as radiation, we have 
investigated the SED of the supercritical accretion flow 
with special attention on the photon-trapping effects. The 
present results are summarized as follows. 



(1) For L<2Le (m <30), at which the photon trapping 
does not occur, the peak frequency of the resultant SED 
becomes, at maximum, three times higher than that ex- 
pected by the standard-disk theory. This is a result of 
the transmission of the high-energy photons from deep in- 
side because of smaller opacity for higher-energy radiation. 
This result might resolve the so-called 'too hot accretion 
disk problem' (Makishima et al. 2000). 

(2) For L>2Le (m>30), the photon-trapping effects 
are substantial. Then, the high-energy photons are more 
efficiently trapped in the flow than the low-energy ones, 
since the former arise near the equatorial plane. Thus, the 
SED shifts towards the lower-energy side as m increases. 

(3) By comparing with the observational data of ULXs 
in the X-ray H-R diagram, we found that the photon- 
trapping effects would only moderately appear in the ac- 
cretion flow of ULXs. We can explain the observed time 
variations of IC 342 SI on the X-ray H-R diagram in terms 
of m modulation. Its black-hole mass is estimated to be 
M ~ 100M©, which is greater than the previous esti- 
mate based on the slim-disk model, M ~ 3OM (Watarai, 
Mizuno, & Mineshige 2001). 

The calculations were carried out at Yukawa Institute 
for Theoretical Physics in Kyoto University. This work is 
supported in part by Research Fellowships of the Japan 
Society for the Promotion of Science for Young Scientists, 
02796 (KO) and 01680 (KW), and the Grants-in Aid of 
the Ministry of Education, Science, Culture, and Sport, 
13640238 (SM). 
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